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Abstract -We employ the random matrix theory (RMT) framework to revisit the distribution 
of resonance widths in quantum chaotic systems weakly coupled to the continuum via a finite 
number M of open channels. In contrast to the standard first-order perturbation theory treatment 
we do not a priory assume the resonance widths being small compared to the mean level spacing. 
We show that to the leading order in weak coupling the perturbative distribution of the 
resonance widths (in particular, the Porter-Thomas distribution at M = 1) should be corrected 
by a factor related to a certain average of the ratio of square roots of the characteristic polynomial 
(‘spectral determinant’) of the underlying RMT Hamiltonian. A simple single-channel expression 
is obtained that properly approximates the width distribution also at large resonance overlap, 
where the Porter-Thomas result is no longer applicable. 


Introduction. Scattering of both classical and quan¬ 
tum waves in systems with chaotic intrinsic dynamics 
is characterised by universal statistical properties m- 
Those can be understood by studying the properties of 
quasi-stationary states in an open system formed at the 
intermediate stage of the scattering process 0H7] . Such 
states are associated with the resonances which corre¬ 
spond to the complex poles £ n = E n — ^T n of the S 
matrix. The resonance widths T n > 0 bring fundamen¬ 
tally new features to the open system in comparison with 
its closed counterpart. In particular, width fluctuations 
govern the decay law ElEO and give rise to nonorthogonal 
modes [HHCLI] leading to enhanced sensitivity to pertur¬ 
bations [12]. Various aspects of lifetime and width statis¬ 
tics are a subject of intensive research, both theoretically 
and experimentally, with recent studies motivated by prac¬ 
tical applications including optical microresonators na. 
superconductor superlattices [14] , many-body fermionic 
systems m, microwave billiards HUE], and dissipative 
quantum maps [I81ll9j . as well as by long-standing inter¬ 
est in superradiance-like “resonance trapping” phenomena 
mmmm- Recent advances in experimental techniques 
have made it possible to test many of theoretical predic¬ 
tions with unprecedented accuracy [2311301 , giving further 


impetus for the theory development. 

The properties of resonance states depend crucially on 
the way in which the system is coupled to the contin¬ 
uum. In the case of the big (“classical”) opening, the wave 
escape follows ballistic trajectories and is dominated by 
system-specific characteristics of chaotic dynamics. The 
quasibound states then have the support in the associ¬ 
ated classical repcller [31], i.e. the (fractal) set of rays 
that remain trapped in the scattering region, and thus 
live long enough to develop random-wave properties 132 . 
The cloud of resonances is formed such that it is typically 
separated by a finite gap from the real axis, with the gap 
width being of the order of the classical decay rate [33] (see 
also [341135] ). Finding resonances inside that gap is a rare 
event, see [19] for a recent review and further discussion. 
Inside the cloud the resonance distribution mainly follows 
the so-called fractal Weyl law Em!, which is determined 
by the fractal dimension of the chaotic repcller. This and 
related characteristics have attracted considerable atten¬ 
tion recently both in physics iT9ll36j and mathematical 
literature El- 

In many realistic situations, however, the opening is 
small (“quantum”) in the sense that it does not scale with 
the system size and supports at a given energy only a finite 
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number M of scattering channels. Representative exam¬ 
ples are microwave cavities [lj (where M is just the number 
of the attached antennas) and quantum dots [38] (then M 
is fixed by transverse quantisation for the modes propa¬ 
gating in the leads). Under such a condition, the waves re¬ 
main in the system long enough to fully explore the whole 
available phase-space and to develop strong diffraction and 
interference effects, leading to decay rates r n being of the 
order of the mean level spacing A. As a result, the scatter¬ 
ing matrix and related observables exhibit statistical fluc¬ 
tuations on the scale ~ A, which are universal, i.e. system- 
independent, provided that the appropriate natural units 
have been used [39]. The only control parameters are then 
the so-called transmission coefficients T c = 1 — |5 CC | 2 , with 
S being the average ‘optical’ S matrix. They quantify the 
system openness, so that T c = 1 (T c 0) stands for a 
perfectly open (almost closed) channel c = 1,..., M. 

The case of the weakly open system, all T c -C 1, is 
of particular interest, since it is most easily realised ex¬ 
perimentally. The resonances are then well separated, so 
one can use perturbation theory 23 to express the widths 
r™ = £ c IK?I 2 < A in terms of the partial amplitudes V£ 
to decay into channel c. The latter have Gaussian statis¬ 
tics inherited from chaotic wavefunctions of the closed sys¬ 
tem. For equivalent channels, the distribution of the res¬ 
onance widths re n = F n /r measured in units of the mean 
partial width (jV/’l 2 ) = F is readily found to be given by 
the x 2 distribution with M/3 degrees of freedom, 




(/?/2) M/3/2 M0/ 2-1 ~/W2 
T(M/3/2) 


(1) 


where /3 = 1 (/3 = 2) indicates the systems with preserved 
(broken) time-reversal symmetry. Distribution © has the 
mean value (re) = M and variance var(re) = (re) , 
showing that typical resonances are indeed very narrow for 
F/ A 1. The regime of overlapping resonances is clearly 
beyond the reach of ©• In particular, it fails to comply 
with the exact Moldauer-Simonius sum rule [411142] for the 
average width (r) Mg = — — T c ). Note that it 
has a log divergence as T c —> 1 related to the emerging 
universal power law far tail ~ T” 2 :43] of the width distri¬ 
bution in such a regime. Describing a crossover to such 
behavior requires going beyond the perturbative scheme 
used to derive © and will be discussed below. 

For a single open channel, distribution © at [3 = 1 was 
first obtained in the classical work by Porter and Thomas 
[40] and bears their name (PTD). It was found to agree 
favorably with early experiments on elastic scattering of 
slow neutrons off heavy nuclei. In fact, for a long time 
that was considered among most firmly established man¬ 
ifestations of chaos in compound nuclear scattering [7]. 
However, recent experiments with improved accuracy [25] 
found evidence of significant deviations from the PTD. 
Various physical mechanisms possibly contributing to the 
discrepancy were discussed in the literature [15I44IH5] - and 
motivated a few proposals towards its resolution. In par¬ 
ticular, one of the proposals [46] suggested a semi-heuristic 


way of going beyond the simple first-order perturbation 
theory for M = 1 and /3 = 1, in order to get access to 
the improved description of those (rare) resonances whose 
reduced widths r„/A ~ 1 despite the weak coupling limit. 

The main goal of our letter is to demonstrate that the 
existing body of knowledge already available (though not 
explicitly worked out) in the literature about the (3 = 1 
case allows one to derive the required expression in a fully 
controllable way. For M = 1 we obtain a simple explicit 
expression ©! that properly approximates the width dis¬ 
tribution even at r n A. Comparing our result to the 
heuristic proposal |46| is instructive, and reveals some in¬ 
teresting similarities and important differences. We also 
reconsider this case from a perspective which is close (but 
not identical) to the spirit of [46] . This allows us to give 
a RMT meaning to the correction factor that arises. 

RMT formulation. - As is well known mu (see also 
[3147] for recent reviews), the non-perturbative approach 
to address universal statistics of N resonance states cou¬ 
pled to M channels is based on applying RMT to the ef¬ 
fective non-Hermitian Hamiltonian of the open system 


U-iV =H- ■ ( 2 ) 

Here H is the Hermitian NxN matrix that corresponds 
to the Hamiltonian of the closed quantum chaotic system, 
whereas the rectangular NxM matrix V consists of the 
coupling amplitudes between N energy levels of H and 
M open scattering channels. The poles of the S matrix 
are just given by N complex eigenvalues of T-LcS which 
therefore coincide with the resonances £ n . Their universal 
fluctuations occur in the RMT limit N 1 (but M kept 
finite) at the local scale A ~ jj, which can be modelled 
by random H drawn from Gaussian orthogonal (GOE, 
f3 = 1) or unitary (GUE, (3 = 2) ensemble, depending 
on the symmetry under time-reversal. Without loss of 
generality, one can restrict the consideration to the spec¬ 
trum center, where the mean level spacing A = Xtt/N 
and 2A is the semicircle radius (A needs to be rescaled if 
E ^ 0). As long as M <C A, the results do not depend 
on statistical assumptions on the coupling amplitudes [35] , 
which are usually taken either as M random [4] or fixed 
[48] A-component vectors characterised by the auxiliary 
coupling constants y c = 7r|| V C || 2 /(2AA). The latter enter 
final expressions only through the transmission coefficients 
T c = 47 c /(1 -I- 7 C ) 2 < 1, so that y c < 1 (y c = 1) corre¬ 
sponds to weak (perfect) coupling. It is also convenient to 
use the renormalized coupling strengths defined as 


9 c = 




> 1 . 


(3) 


The standard treatment of the regime of weak cou¬ 
pling to the continuum amounts to considering the anti- 
Hermitian part of "Heff as a perturbation of the Hermitian 
part H. To the leading order in y c -C 1, the resonance 
widths are given by r„ = (n\VV^\n) = X^cl^nl 2 ’ w h ere 
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| n) stands for the eigenvector corresponding to the n-th 
eigenvalue (energy level) of H. A rotation that diagonal¬ 
izes random H ensures that in the limit N 1 the compo¬ 
nents become independent Gaussian distributed vari¬ 
ables. For equivalent channels this readily gives ©, with 
F = 27 A/- 7 T « TA/27T. (The generalisation for nonequiv¬ 
alent and correlated weak channels is also known [15]. 1 
A non-pertubative approach to the resonance distribu¬ 
tion is based on making use of an electrostatic analogue 
[HI51I35] treating the S matrix poles as two-dimensional 
charges and employing the powerful supersymmetry tech¬ 
nique [IQ to perform the statistical average. Let us stress 
that in such a way the distribution of the resonance widths 
was already derived for any number of arbitrary open 
channels first at /3 = 2 mm and then at B = 1 [50] - Un¬ 
fortunately, in such generality the exact /3 = 1 expression 
takes a rather complicated form of a three-fold integral and 
is not very transparent, calling for further analysis. Intro¬ 
ducing the coupling strengths g c defined in ([3]) and arrang¬ 
ing them in the increasing order g = g\ < c /2 < • ■ • < <?m ) 
the exact probability density for the distribution of scaled 
widths y n = 7rr n /A > 0 reads as follows [50] : 

d 2 

Pm(v) = -j~2 Pm(v ), (4) 

where the ‘potential’ function Fm{u) is given by 


Pm(u) 


I 

7r 


f°° dp 2 X2{.P2)e z p2 


Im(P2,V ) (5) 


/m 


f P2 d Pl (p 2 - 

■A y/Pi - 


<t>M(jpi,P2,y) 



Pi) Xi(Pi) e 2Pl 

"i U^iVW^l 


<t>M{pi,P2,y) 


dA ri _ eVX ltLi(9c-X) 

2 ’(\-Pin\-p2) 2 ’ 


( 6 ) 

(7) 


with the factors xi(p) and X 2 (p) each showing a rather 
complicated alternating pattern of values ±1 or 0 for var¬ 
ious values of p, see [50] . For the goals of our discussion 
it is enough to know only that \i (p) = 1 for 1 < P < g, 
whereas \2 (p) = 0 in the same range. This distribution 
satisfies the Moldauer-Simonius relation exactly, yielding 



c= 1 


f g c + 1 

\9c - 1 


( 8 ) 


for the average width. Note also that the validity of the ex¬ 
act distribution (]4]) at M = 1 was verified in high-precision 
scattering experiments with microwave billiards [23] . 

Weak coupling approximation. - Our goal now is to 
investigate the weak coupling limit g = g\ ^$> 1 of the exact 
expression considering, however, arbitrary values of the 
widths y. When performing such a limit we use that due 
to the nature of \i, 2 (p) factors the leading contribution to 
the integrals comes from the domain (A,pi) ~ 1 Cg <P 2 - 
Correspondingly, we can replace factors (A — P 2) 2 ~ p\. 


P 2 — Pi ~ P 2 , and \/p\ — 1 ~ P 2 as well as approximate 
g c — pi « g c and g c — A ~ g c for all c = 1,... M. Further 
analysis also shows that one can replace the upper limit 
P 2 in the integral d 6 | by infinity. In this way we arrive at 


Pm{v) 


where 



dp 2 X 2 {P 2 ) e ? P2 

pi n^ 1 ii-5c 1 p 2 i i /2’ 


Hy) 



dpi e 2 pi 
v^T(A-pi ) 2 


(9) 


( 10 ) 


After somewhat lengthy calculations one is able to evalu¬ 
ate the double integral (flOl) in a closed form, finding 


$(y) = 2 L 


-ft'o(l) (coshy - bmh;/ ) + R'i(l) sinh y 


( 11 ) 


where K 0: i{x) are the modified Bessel functions. 

Noticing that 4>(y) in © is a ^-independent factor and 
taking into account also that P 2 > g due to X '2 ip 2 ) function 
in the integrand, one may naively conclude that for g > 1 
the leading order approximation to the width distribution 
in Q is obtained by twice differentiating the g-dependent 
integral factor in ([5]). This results in 


P { M aive \y) = Hy)V x 2 (y), 


where 


p x 2 (v) 


f°° dp \2 ip) e ^ 

A ^n^iii-sc-vi 1 / 2 ’ 


( 12 ) 


(13) 


Note that expression m has a specific meaning: it is 
just a version of the probability density of the \ 2 distribu¬ 
tion for weak non-equivalent channels 1491150] . which can 
be shown to reduce to the standard form m for all the 
same g c = < 7 , with g « A. a t 7 -C 1. Therefore, all the 
nontrivial leading corrections to the PTD in such an ap¬ 
proximation seem to be accumulated in the factor <h(y). 
However the following problem arises: keeping only such 
leading term results in losing the correct normalization, 
being thus an inconsistent approximation to the probabil¬ 
ity density. Thus the other (higher-order) terms need to 
be also considered, making the analysis more involved. 

A resolution to that problem can be found by analysing 
and exploiting the factorised structure of j9| further. That 
form suggests an Ansatz providing the leading-order ap¬ 
proximation to the width distribution (|4]) by differentiat¬ 
ing the g-dependent factor only once as follows: 


Pm{p) = - A [$(y)W M (y)] . (14) 

Function Am (y) has the meaning of the cumulative (tail) 
distribution function associated with the density (fl3l) . be¬ 
ing normalized as Am( 0 ) = 1 and explicitly given by 


Am( 2 /) 


r°° dp X 2 (p)e 2 P 

A ^n"iii -9c-v /2 ' 


(15) 
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Fig. 1: Distribution of the scaled resonance widths (in units of 
the mean width) for single-channel scattering at g = 1.1 (•), 
1.2(o), 1.5 (★), and 3(D). The solid lines stand for approxi¬ 
mation m whereas the dashed line is the Porter-Thomas ex¬ 
pression (M = 1). The symbols correspond to numerics with 
10 5 realizations of 250 x 250 GOE random matrices (only 25 
levels around the spectrum center were taken for each real¬ 
ization). The inset shows the comparison of the mean width 
(y) obtained with Eq. (1161) to the exact (Moldauer-Simonius) 
expression ( y) MS as a function of the coupling strength g. 


Note that the form (fTdl) satisfies the normalisation condi¬ 
tion automatically. To ensure a positive probability den¬ 
sity, the function 4>(2 /)A/m( 2/) needs however to be mono¬ 
tonic. We can show that this is indeed the case for, e.g., 
M = 1,2 where our Ansatz is therefore fully consistent 
In particular, at M = 1 we arrive after the straightforward 
integration to an attractively simple formula 


^ > M=i(y) 


d 

dy 


$(j/) erfc 



(16) 


where erfc(z) = 1 — e t dt is the complementary 

error function. The formulae (fill) , (fTdl) (fl6l) represent the 
consistent leading-order approximation of the resonance 
width distribution at weak coupling for any value of y : 
and as such constitute the main result of this paper. 

It is instructive to analyse the single-channel case in 
more details. In the region of small widths, j/ C 1, dis¬ 
tribution (fl6l) coincides with the PTD up to the terms of 
the order ~ 0(y 3 / 2 ). It shows, however, significant devi¬ 
ations for large widths, with its leading asymptotic being 

Vi(y » 1) = e-to-W^iT^ ^^ ^+Ofe- 3 )]. 

One sees that by formally setting g = 1 here, we get the 
power-law behavior ss l/i/ 2 y 2 , which up to a constant 
reproduces the known exact one « 1/2 y 2 . Correspond¬ 
ingly, the average width ( y ) computed by using (1161) de- 

1 For M > 4 the positivity of the ensuing density turns out to 
be mildly violated for a relatively small region of narrow widths. It 
can possibly be interpreted as a precursor of the depletion of narrow 
resonances due to the gap formation at larger M, see EH- Formally, 
it indicates on the noncommutative limits M 1 and g 1, thus 
nonuniform convergence due to the approximations involved. 


velops ^=| ln (<7 — 1 )| divergence at g —> 1 , in qualitative 
agreement with the exact expression (y ) Ms given by ©. 
This suggests that even though our formula (fl 6 l) has been 
derived under the assumption of g 1 , it should prop¬ 
erly approximate the exact distribution also at moderate 
g ~ 1 , provided that the rescaled widths y/ (y) are used. 
To verify this procedure, we have performed numerical 
simulations with the GOE matrices and indeed found a 
good agreement, as is clearly demonstrated in fig. 1 . 

We briefly mention that similar findings hold at M = 2 
open channels (not shown here), function (fl5l) being then 
given by Af 2 (y) = y/gTgi &[(g 2 ~p)(p-gi)}~ 1/2 e~% p . 

Relation to spectral determinants. To shed fur¬ 
ther light on the meaning and origin of the correction fac¬ 
tor ^(y), we note that it is independent of M and g and 
thus can be studied for the particular case M = 1. We 
now reconsider this case by a different method inspired, 
in part, by our attempts to understand the nature of the 
heuristic proposal [56]. We start similarly with the ex¬ 
act joint probability density V{£ k } of N resonances in the 
complex plane, but find it slightly more convenient to work 
with the ‘fixed coupling’ version of the model EH: 


N 


V{£ k } oc 7 n S(£ n ) n |g r g : 12 <5 


m<n 


| £m ~ £n 



(17) 

where S(£ n ) = ( 7 ^) 3 / 2 exp[-^( 7 2 A 2 + £ 2 - T 2 /4)] 
and we have omitted numerical constants which may de¬ 
pend on N but are otherwise parameter-independent. 
Making use of this expression and following the method 
proposed in m for the /3 = 2 case of the model (see [55] 
for a detailed consideration), one can arrive in the limit 
N ^$> 1 to a convenient representation for the ensuing den¬ 
sity p(E, T) of resonances around the point £ = E — |I’. 
Note that at the spectrum center, E = 0, this function 
has precisely the meaning of the corresponding resonance 
widths distribution. Under a weak condition T < 1 but 
otherwise for arbitrary resonance overlap , we find 


Pu{v) =is. v K 0, 77) “ 71 <r "' 2 ' M!,) ■ (18) 


where v k 7 ( 2 /) is the following spectral determinant ratio: 


*7 (y) 


lim Cn 

N—^oo 


|det(^- 1) +z2/A/2V)| 2 \ 

I det("Hgg _1)t + iy\/2N)\ J 


(19) 


Here the angular brackets stand for the GOE average over 
the Hermitian part of 1 , which has now the reduced 
dimension of N — 1. The constant Cn is indicated for en¬ 
suring the limit is well defined. We stress that expressions 
(E| El are the exact representation of the resonance 
density at M = 1 in the limit considered, being valid for 
any values of width y and at arbitrary coupling strength 
g > 1. The above representation is especially convenient 
for our purposes, since it factorizes into a product of the 
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PTD factor -J= e sy / 2 and the factor 4/ 7 (y), which there¬ 
fore accumulates all corrections beyond the PTD. 

Unfortunately, performing exactly the GOE average in 
ED at arbitrary coupling 7 ~ 1 is a rather challenging 
task and remains an interesting outstanding problem of 
RMT. However, in the weak-coupling limit 7 C 1 (i.e. 
j > 1 ) we can easily find a consistent approximation 
to that factor at the leading order in 7. This obviously 
amounts to simply setting 7 = 0 inside the RMT average, 
that is to neglecting the anti-Hermitian part by replacing 
^ with (N — 1) x ( N — 1) matrix H from GOE in 
m- After doing this, we arrive at the problem of evalu¬ 
ating the following GOE average 

lim $~(|/)= lim C N ( ^det (H 2 + y 2 A 2 /4iV 2 )\ , (20) 

where now Cn = 77 ^ 7 ^- Actually, the right-hand side 
in (l2l)l) is similar to the objects appeared already in other 
RMT contexts, being the subject of a separate work [551 . 
A particularly pleasant though not unexpected outcome of 
that work is that (1201) can be evaluated in a closed form, 
yielding lim 7 = o'k 7 (j/) = <f>(y) where <h(y) is given exactly 
by ED. This independently confirms the calculation done 
in the first part of this paper for M = 1 by a completely 
different method and thus the internal consistency of all 
approximations made in the course of our derivation. 

Discussion and conclusions. - We are now in the po¬ 
sition to make contact to the proposal of ref. [46]. Basi¬ 
cally, that proposal coincides with the form of the “naive” 
(not-normalized) expression (1121) in which the correct fac¬ 
tor $( 2 /) was however replaced by a heuristic expression 
( sm ^ v ) ■ Although such a form at (3 = 1 shares sim¬ 
ilar larg e-y asymptotic with $(?/), it is manifestly differ¬ 
ent from the latter, rendering the proposal [55], in the 
strict sense, invalid. At this point it is appropriate to 
mention that in [55] the authors actually arrived, though 
in a somewhat different way, to the problem of evalu¬ 
ating the same average featuring in the right-hand side 
of (HOD . However, they proposed to replace true GOE 
spectra there by a picket-fence model of equidistant lev¬ 
els. Though one may hope that the phenomenon of level 
repulsion could make it qualitatively sensible, such a re¬ 
placement is certainly not a systematic approximation. In 
particular, neglecting spectral fluctuations is known to un¬ 
derestimate interference effects related to non-orthogonal 
wavefunctions in open systems [551I551 . It is interesting to 
note, however, that similar manipulations with the exact 
GUE formulae for the resonance width distribution [43j 
show that the “true” correction factor is actually given by 
limjv^oo (det (Df 2 + y 2 X 2 /4N 2 )) oc S 111 J 1 y , thus coinciding 
with the proposal [46] at /3 = 2. Such a coincidence ap¬ 
pears nevertheless as a stroke of luck, though somewhat 
stronger level repulsion in the GUE case may indeed con¬ 
tribute to its success. 

In conclusion, we used RMT framework to develop a 
systematic approximation to the resonance width distri¬ 


bution in weakly open quantum chaotic systems beyond 
the standard first-order perturbation theory. For the small 
number M = 1,2 of open channels we found simple ex¬ 
pressions which accurately reproduce the width distribu¬ 
tion in the regime of moderately overlapping resonances, 
remaining qualitatively correct even at strong coupling to 
the continuum when resonances overlap considerably. As 
such, our results may be advantageous for practical appli¬ 
cations. We also believe that the new representation (H51) 
may prove to be useful for finding alternative derivation 
of the exact distribution in a simpler explicit form. Com¬ 
parison with similar expressions at /? = 2 in [52] makes us 
to anticipate that the factorized form (1181) may retain its 
validity for M > 1 after replacing the PTD factor with 
the density (fl3l) . It remains an outstanding challenge to 
verify this from the first principles. Finally, it is worth 
mentioning that there is growing interest in this type of 
problems also in mathematical literature, see e.g. .551. 

* * * 

We are grateful to V.V. Sokolov for numerous dis¬ 
cussions. YVF acknowledges support by EPSRC grant 
EP/J002763/1 “Insights into Disordered Landscapes via 
Random Matrix Theory and Statistical Mechanics”. 

REFERENCES 

[1] Stockmann H.-J., Quantum Chaos: An Introduction 
(Cambridge University Press, Cambridge, UK) 1999. 

[2] Kuhl U., Stockmann H.-J. and Weaver R., J. Phys. 
A: Math. Gen., 38 (2005) 10433. 

[3] Gradoni G., Yeh J.-H., Xiao B., Antonsen T. M., 
Anlage S. M. and Ott E., Wave Motion, 51 (2014) 606. 

[4] Sokolov V. V. and Zelevinsky V. G., Nucl. Phys. A, 
504 (1989) 562. 

[5] Fyodorov Y. V. and Sommers H.-J., J. Math. Phys., 
38 (1997) 1918. 

[6] Fyodorov Y. V., Savin D. V. and Sommers H.-J., J. 
Phys. A: Math. Gen., 38 (2005) 10731. 

[7] Mitchell G. E., Richter A. and Weidenmuller 
H. A., Rev. Mod. Phys., 82 (2010) 2845. 

[8] Savin D. V. and Sokolov V. V., Phys. Rev. E, 56 (1997) 
R4911. 

[9] Dittes F.-M., Phys. Rep., 339 (2000) 215. 

[10] Schomerus H., van Bemmel K. J. H. and Beenakker 
C. W. J., Europhys. Lett., 52 (2000) 518. 

[11] Savin D. V., Legrand O. and Mortessagne F., Euro¬ 
phys. Lett., 76 (2006) 774. 

[12] Fyodorov Y. V. and Savin D. V., Phys. Rev. Lett., 108 
(2012) 184101. 

[13] Schomerus H., Wiersig J. and Main J., Phys. Rev. A, 
79 (2009) 053806. 

[14] Gluck M., Kolovsky A. R. and Korsch H. J., Phys. 
Rep., 366 (2002) 103. 

[15] Celardo G. L., Auerbach N., Izrailev F. M. and 
Zelevinsky V. G., Phys. Rev. Lett., 106 (2011) 042501. 

[16] Poli C., Legrand O. and Mortessagne F., Phys. Rev. 
E, 82 (2010) 055201. 


P-5 








Yan V. Fyodorov 1 Dmitry V. Savin 2,3 


[17] Weich T., Barkhofen S., Kuhl U., Poli C. and 
Schomerus H., New J, 16 (2014) 033029. 

[18] Keating J. P., Novaes M., Prado S. D. and Sieber 
M., Phys. Rev. Lett., 97 (2006) 150406. 

[19] Novaes M., J. Phys. A: Math. Theor., 46 (2013) 143001. 

[20] Persson E., Rotter I., Stockmann H.-J. and Barti-i 
M., Phys. Rev. Lett., 85 (2000) 2478. 

[21] Auerbach N. and Zelevinsky V., Rep. Prog. Phys., 74 
(2011) 106301. 

[22] Mulhall D., Phys. Rev. C, 91 (2015) 014305. 

[23] Kuhl U., Hohmann R., Main J. and Stockmann H.- 
J., Phys. Rev. Lett., 100 (2008) 254101. 

[24] Dietz B., Friedrich T., Harney H. L., Miski-Oglu 
M., Richter A., Schafer F. and Weidenmuller 
H. A., Phys. Rev. E, 81 (2010) 036205. 

[25] Koehler P. E., Becvar F., Krticka M., Harvey 
J. A. and Guber K. H., Phys. Rev. Lett., 105 (2010) 
072502. 

[26] Di Falco A., Krauss T. F. and Fratalocchi A., Appl. 
Phys. Lett., 100 (2012) 184101. 

[27] Liu C., Di Falco A. and Fratalocchi A., Phys. Rev. 
X, 4 (2014) 021048. 

[28] Barkhofen S., Weich T., Potzuweit A., Stockmann 
H.-J., Kuhl U. and Zworski M., Phys. Rev. Lett., 110 
(2013) 164102. 

[29] Kuhl U., Legrand O. and Mortessagne F., Fortschr. 
Phys., 61 (2013) 404. 

[30] Gros J.-B., Kuhl U., Legrand O., Mortessagne F., 
Richalot E. and Savin, D. V., Phys. Rev. Lett., 113 
(2014) 224101. 

[31] Lu W. T., Sridhar S. and Zworski M., Phys. Rev. 
Lett., 91 (2003) 154101. 

[32] Schomerus H. and Tworzydlo J., Phys. Rev. Lett., 93 
(2004) 154102. 

[33] Gaspard P. and Rice S. A., J. Chem. Phys., 90 (1989) 
2225. 

[34] Haake F., Izrailev F. M., Lehmann N., Saher D. 
and Sommers H.-J., Z. Phys. B, 88 (1992) 359. 

[35] Lehmann N., Saher D., Sokolov V. V. and Sommers 
H.-J., Nucl. Phys. A, 582 (1995) 223. 

[36] Altmann E. G., Portela J. S. E. and Tel T., Rev. 
Mod. Phys., 85 (2013) 869. 

[37] Nonnenmacher S., Sjostrand J. and Zworski M., 
Ann. Math., 179 (2014) 179. 

[38] Alhassid Y., Rev. Mod. Phys., 72 (2000) 895. 

[39] Guhr T., Muller-Groeling A. and Weidenmuller 
H. A., Phys. Rep., 299 (1998) 189. 

[40] Porter C. E. and Thomas R. G., Phys. Rev., 104 
(1956) 483. 

[41] Moldauer P. A., Phys. Rev., 171 (1968) 1164. 

[42] Simonius M., Phys. Lett. B, 52 (1974) 279. 

[43] Fyodorov Y. V. and Sommers H.-J., JETP Lett., 63 
(1996) 1026. 

[44] Weidenmuller H. A., Phys. Rev. Lett., 105 (2010) 
232501. 

[45] Volya A., Phys. Rev. C, 83 (2011) 044312. 

[46] Shchedrin G. and Zelevinsky V., Phys. Rev. C, 86 
(2012) 044602. 

[47] Fyodorov Y. V. and Savin D. V., Resonance scattering 
of waves in chaotic systems , in The Oxford Handbook of 
Random Matrix Theory, edited by Akemann G., Baik 
J. and Di Francesco P., (Oxford University Press, UK) 


2011, pp. 703-722 arXiv: 1003.0702 . 

[48] Verbaarschot J. J. M., Weidenmuller H. A. and 
Zirnbauer M. R., Phys. Rep., 129 (1985) 367. 

[49] Alhassid Y. and Lewenkopf C. H., Phys. Rev. Lett., 
75 (1995) 3922. 

[50] Sommers H.-J., Fyodorov Y. V. and Titov M., J. 
Phys. A: Math. Gen., 32 (1999) L77. 

[51] Fyodorov Y. V. and Khoruzhenko B. A., Phys. Rev. 
Lett., 83 (1999) 65. 

[52] Fyodorov Y. V. and Sommers H.-J., J. Phys. A: Math. 
Gen., 36 (2003) 3303. 

[53] Fyodorov Y. V. and Nock A., J. Stat. Phys., 159 
(2015) 731. 

[54] Poli C., Savin D. V., Legrand O. and Mortessagne 
F., Phys. Rev. E, 80 (2009) 046203. 

[55] Savin D. V. and De Vaulx J.-B., Acta Phys. Pol. A, 
124 (2013) 1074. 

[56] Killip R. and Kozhan R., arXiv:1501.05160 (2015). 


P-6 



